Draft version February 2, 2008 

Preprint typeset using I^T^jX style cmulatcapj v. 6/22/04 



O 
O 

(N 



> 

(N 
t> 
O 

o 

Or 
i 

o 
C3 



AN AB INITIO APPROACH TO SOLAR CORONAL LOOPS 

B.V. GUDIKSEN AND A. NORDLUND 
Draft version February 2, 2008 

ABSTRACT 

Data from recent numerical simulations of the solar corona and transition region are analysed and the 
magnetic field connection between the low corona and the photosphere is found to be close to that of 
a potential field. The fieldline to fieldline displacements follow a power law distribution with typical 
displacements of just a few Mm. Three loops visible in emulated Transition Region And Coronal 
Explorer (TRACE) filters are analysed in detail and found to have significantly different heating rates 
and distributions thereof, one of them showing a small scale heating event. The dynamical structure 
is complicated even though all the loops are visible in a single filter along most of their lengths. None 
of the loops are static, but are in the process of evolving into loops with very different characteristics. 
Differential Emission Measure (DEM) curves along one of the loops illustrate that DEM curves have 
to be treated carefully if physical characteristics are to be extracted. 

Subject headings: Sun: corona - Sun: magnetic fields - MHD 



1. INTRODUCTION 

Since the launch of the SOlar and Heliospheric Ob- 
servatory (SOHO) and later the Transition Region And 
Coronal Explorer (TRACE), observations of EUV loops 
in the solar corona have become sufficiently good to 
provide a te s tbed for coronal loop models. Earlier, 
iRosner et alJ l|1978j) produced scaling relations for loops 
with constant pressure, uniform heating and constant 
cross section along their lengths, the so calle d RTV scal- 
ing la ws. These were later generalized by iSerio et"al~l 
Ij 198 1(1 . including non-uniform heating, and non-uniform 
pressure in the form of two scale heights. It has though 
become increasingly clear that even these hydrostatic 
models have problems reprod ucing the majority of the 
loops observed in the EUV. lAschwanden et all 1)20011) 
compared 41 loops observed with TRACE with stan- 
dard loop models, and found that none of the observed 
loops could be fitted with the RTV scaling laws and only 
roughly 30 % of the loops could be explained by hydro- 
static soluti ons with foot point heatin g. 

Recently iWinebareer et aD l)2003f) compared both 
EUV and X-ray loops and found that long EUV loops are 
overdense by up to 3.4 orders of magnitude compared to 
hydrostatic uniformly heated models. Furthermore, only 
28 % of the 67 loops were explainable by hydrostatic non- 
uniformly heated models. This included X-ray loops with 
an assumed filling factor of unity, which could in principle 
be reconciled with the models by assuming a smaller fill- 
ing factor, but this in itself would introduce other prob- 
lems. Reducing and interpreting EUV coronal loop ob- 
servations are complicated and cannot provide all the 
necessary information. Obtaining reliable estimates of 
gas parameters and velocities for a whole loop demands 
spectral information only available by rastering, which 
compromises the understanding of the dynamic nature 
of loops. 

Further problems concerning proper background sub- 
traction for both spectra and imaging ha ve been shown 
to have effects on the ded u ctions made llMarte ns et al.l 
120021 ISchmelz et all 120031 iDel Zanna fc Masonl 12003ft . 
Modelling of observed loops consequently incorporates 



assumptions not necessarily confirmed by observations. 
Common for most models is the assumption of a uniform 
loop cross section. This assumption is unchallenged by 
observations since the measured expansion of the loops 
is severely limited by the spatial resolution of existing in- 
struments. Lately, claims have bee n made that loops now 
appear to be resolved by TRACE <|Testa et al.ll2002ft . At 
the same time, loops of cross section on the order of a 
few pixels show signs of consisting of smaller structures 
since estimates of cooling time, based on the evolution 
of the intensity in two TRAC E pass bands, do not agree 
with hydrostatic loop models ijWarren et al.ll200^) . Re- 
sults from a 3D MHD model of the corona gives us the 
opportunity to compare individual loops identified in the 
model with observations and investigate if the often used 
approximations in loop models are valid. 

The problems of modelling coronal loops will here 
be treated with data f r om th e numerical simulations of 
Gudi ksen fc Nordlundl l)2004|) . We will discuss the mag- 
netic field state and the validity of a potential field ex- 
trapolation of a magnetogram, the assumption of con- 
stant circular cross sections for loops, as well as looking 
at flows, energy balance and gas parameters for a few 
single loops. The time evolution of one of these loops 
and the DEM curves for a few points along the same 
loop will also be discussed. Emission, Doppler shifts, 
and non-thermal widths of a number of selected spectra l 
lines are treated in a separate paper l|Peter et alJ l2004). 

2. MODEL DESCRIPTION 

The numerical simulation of the solar corona 
is described in detail in a separate paper 
ijCudiksen fc Nordlundl I2004D and is only introduced 
briefly here. The model is based on an MHD code 
that incorporates relevant physics, including a radiative 
cooling function and Spitzer conductivity. The model 
spans 60 x 60 Mm 2 of the solar surface, and stretches 
37 Mm up into the corona from the photosphere. A 
model for the solar photosphcric velocity field stresses 
the magnetic field which initially is a potential extrapo- 
lation of an MDI/SOHO high resolution magnetogram 
of AR 9114, scaled to fit inside the box. The initial 
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Fig. 1. — Histogram showing the distribution of horizontal dis- 
tances a field line travels from the photosphere to the lower corona 
for a typical snapshot (full) and a potential extrapolation of the 
same field (dashed) 



thermal stratification is a FAL-C llFontenla et aljH99l 
atmosphere in the photosphere and chromosphere 
with an isothermal 1 MK corona above the transition 
region. The temperature in the lower, optically thick 
atmosphere is forced toward the initial temperature 
profile on a typical timescale of O.f s in the photosphere 
and decreasing as p 2 ^ 3 , making it unimportant in the 
transition region and corona. To check the effect of 
the chromospheric stratification on the corona a "cool" 
chro mosphere with no av erage temperature increase 
(see iCarlsson fc Steinll20f)2l and references therein) was 
also used. The choice of chromospheric model is not of 
great importance for the loop structures, and only the 
model with the standard FAL-C chromosphere will be 
discussed here. 

3. RESULTS 

The definition of a loop is often unclear, and can lead 
to a number of misconceptions, so we will explain what 
we mean by the term "loop" and what the consequences 
are. A loop observed in a narrow wavelength band is de- 
fined by the plasma emitting in that certain wavelength 
band. In principle that definition is independent of the 
magnetic field, but in the solar corona it can, to a very 
good approximation, be assumed that the near isother- 
mal plasma seen emitting is caught in a flux bundle, be- 
cause of the efficient thermal conduction along the mag- 
netic field. Often magnetic field lines are used to trace 
the magnetic field, but it is important to emphasize that 
a field line is not a physical entity, it is simply a line 
following the direction of the magnetic field. It is there- 
fore not affected by changes in the local magnetic field 
strength, and cannot be followed in time. The magnetic 
field in the solar corona is space filling and it is rarely, if 
ever, possible to identify a loop from the magnetic field 
alone, without the information from narrow band imag- 
ing instruments. Figures where loops are shown as a 
number of field lines are made by a subjective choice of 
field lines that along their lengths have gas emitting in 
a particular wave length band. Consequently, loops are 
the result of plasma being near isothermal, and not the 
effect of the magnetic field being special in the volume 
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Fig. 2. — Histogram showing the difference between foot point 
locations for field lines going through the same points in the lower 
corona, for a potential extrapolation and a snapshot 



defined by the emitting plasma. A visible loop is a result 
of the heating history within the flux bundle in which the 
emitting gas is caught and not a result of the magnetic 
field in the flux bundle being very much different from 
the magnetic field nearby. 

It is often assumed that the emitting gas and the same 
magnetic flux are involved in a long lived loop, and even 
though this is often the case, it need not be. In a 3D sim- 
ulation such as this, one is fortunate enough to able to 
choose between identifying loops by following the mag- 
netic field, and identifying them from a high level of emis- 
sion. In case of time evolution following a loop based 
on the magnetic field is not possible — at least not in a 
unique way — so instead one has to follow a plasma parcel 
through time. Only to the extent that the field is "frozen 
in" does this mean that the same flux is passing through 
the plasma parcel. 

3.1. Loop connection from photosphere to corona 

Observed polarisation signals from the Sun usually 
originate from a relatively thi n layer in the photo sphere 
(but see recent results from Solan kTet al J 12003). The 
connection to a 3D magnetic field is then usually made 
by using the observed photospheric magnetogram and ex- 
trapolating by a method producing either a linear force 
free or potential field. With the movements of the mag- 
netic foot points in the solar photosphere that method 
must be questioned. To investigate whether this method 
provides a good approximation we have compared the 
horizontal distance a fieldline traverses from the photo- 
sphere to a predetermined height just above the transi- 
tion region. The magnetic fieldlines are followed from a 
horizontal layer just above the transition region where 
they are distributed evenly over the whole layer. Fig- 
ure Q] is a histogram showing the distribution of the dis- 
tances for a typical snapshot and for a potential extrap- 
olation of the vertical magnetic field in the photosphere. 
Even though the difference between these distributions 
is very small, the distance between the foot point of 
a field line passing through a point in the low corona 
obtained from the simulation and the assumption of a 
potential magnetic field may be very large for individ- 
ual fieldlines. Figure [5] shows the distance between foot 
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Fig. 3. — Foot point widths in x and y at the photosphere (dots) 
and in the transition region (crosses) for loops with a circular cross 
section and a diameter of 0.8 Mm at their top. Foot points of 
bright TRACE loops are encircled 



points of a fieldline passing through the same point in the 
lower corona from the potential extrapolation and from 
the snapshot. Around Qu asi-Separatrix-Layers (QSLs) 
((Priest fe Demouhnl Il995|) the computed distances are 
generally very large, because the shuffling of the fieldline 
foot points make the QSLs move around in the atmo- 
sphere, making fieldlines starting at the same geomet- 
rical point end up at large distances from each other. 
Fieldlines through such points are not included in Fig. [21 
The PDF is consistent with a power law but since we do 
not have the full range of driving scales included, and be- 
cause the simulation has run for only one turnover time 
of the largest driving scale it is plausible that we would 
get a high distance tail if the simulation had run for a 
longer period. For the majority of the fieldlines followed 
here the difference in distance is modest, on the order 
of a few Mm. This is a small change in connectivity, 
and in spite of the maximum driving scale being smaller 
than the super granular scale typically used, substantial 
heating is produced. 

Fieldlines that do not reach above the transition re- 
gion are disturbed the most. This is partly because the 
higher densities in the chromosphere give low propaga- 
tion speeds for disturbances, and thus allow larger distor- 
tions. Correspondingly, roughly 90 % of the dissipated 
energy i s injected in the lower atmosp here and not in the 
corona (iGudiksen fc Nordlundl l2004). Large scale shear 
in the corona can only be created by a large scale per- 
sistent photospheric velocity field, if it is not already a 
property of the magnetic field at the time of emergence. 

3.2. Loop cross section 

Loops have generally been assumed to be of constant 
cross section at all heights in the corona, based on the ob- 
servations from SOHO and TRACE, where there seems 
to be no increase in cross section with height, contrary 
to what would be expected for a potential field. Implic- 
itly, to at all reach conclusions about cross sections from 
such observations, it has been assumed in previous works 
that cross sections of flux surfaces are roughly circular. 
For the heating mechanism evident in this work, such an 
assumption is highly unlikely. To achieve heating, mag- 




FlG. 4.— Emulated TRACE 171 (top) and 195 (middle) with 
an emulated magnetogram (bottom), showing three loops bright 
in emission, and 5 points as crosses in the emulated TRACE 171 
image which are used to create DEM curves 



netic fieldlines have to be at an angle relative to each 
other at the heating location. Identifying a loop by its 
EUV or X-ray emission and therefore indirectly by its 
heating history, one can be certain that the field lines in 
such a loop are not aligned. The cross section of a flux 
surface that is circular at one point along a loop will not 
remain circular in a tangled fieldline geometry. 

Figure shows the horizontal width of loops that have 
a circular-loop cross section at their tops of 0.8 Mm. 
Loops with a wide variety of maximum heights have been 
chosen, but there is no obvious dependence of the changes 
in cross section on the maximum height. It is clear that 
the loops have significantly smaller cross sections in the 
photosphere, and that they are non-circular. When mov- 
ing further up into the transition region, the loops be- 
come wider, some having expanded by a factor of ten 
but generally the expansion is not identical in all direc- 
tions. That makes the traditional "wine-glass" picture of 
the magnet ic field in the transition region a bad approx- 
imation, as lSchriiver fc Titla 1)2003(1 have pointed out af- 
ter modelling of network flux concentrations. There is 
a tendency for bright TRACE loops to be very asym- 
metric in the photosphere. It is clear that if these loops 
are assumed to have circular cross sections at their tops, 
few (if any) have circular cross sections at a height cor- 
responding to the transition region. Therefore neither 
cross section nor cross section area are conserved along 
any of the loops we have followed. In general the foot 



4 



5 



o 

Q_ 
<D 
Q 

>- 
cn 
i_ 

CD 
C 



0.006 r 

0.004 ; 
0.002 7 

o.ooo - 

-0.002 ; 

-0.004 ; 

-0.006 L 

2 



- 15 



» - 1 



1 12 14 16 14 
Z [ Mm ] 



6 4 



- 5 



JO w 
o 

2 < 



Fig. 5. — Loop 1 with left and right leg in Fig. |1] being to the left and right respectively, for the energy deposition showing heating 
(solid), radiative cooling (dashed), Spitzer conductivity (dash-dot-dot-dot), convective energy deposition (long dashes) and emission in the 
TRACE 171 A filter(dotted) as a function of height 
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Fig. 6. — Gas parameters for loop 1 showing density (solid), gas pressure (dashed), temperature (dotted) and the absolute velocity along 
the magnetic field (long dashed). For this loop the gas flows down both legs 



prints of the loops shown in Fig. |31 are wrinkled and not 
geometrically simple, and their area cannot be approxi- 
mated by product of AX and AY. We expect the foot 
prints to increase in complexity with increasing resolu- 
tion. 

It should be pointed out that the loops followed here 
are not chosen on the basis of a visible characteristic, 
but are simply based on choosing random points in the 
corona, whereafter a set of fieldlines around each point 
are followed to the photosphere. This means that not 
all the loops treated here are representative of the ob- 
served cross sections of TRACE loops, but lend evidence 
to the fact that bright loops in general do not have cir- 
cular cross sections along their whole length. In order to 
probe the true cross section of loops, the spatial resolu- 
tion would have to be raised by at least a factor of two, 
to include velocity gradients produced by flow scales in 
the photosphere that the resolution here cannot resolve. 

3.3. Loop Heating 

Loops have in general highly time dependent heating 
profiles, and the heating is independent in the loop legs. 
It is thus only possible to define a general heating pro- 



file if a large number of loops are used. We obtain an 
average heating profile that above the transition region 
decreases exponentially, with a scale length of ~ 6 Mm. 
The height dependence of the heating is highly variable 
from loop to loop. The scale height behaviour emerges 
since the heating is proportional to the current squared 
and because the magnetic field in the corona is close to 
force free, which makes the current along each fieldline 
proportional to the magnetic field. In addition to being 
close to force free, with the property J cx B, the mag- 
netic field also does not deviate much from a potential 
field (cf. Sec. 13.1(1 . The overall magnetic field strength 
thus decreases roughly exponentially with height, with 
a scale height determined by the distance between the 
dominating magnetic polarities. In the present periodic 
model the shortest distance between the two main po- 
larities is roughly half the box size, so at coronal heights 
the wave number corresponding to the box size ko dom- 
inates. Assuming that the magnetic field drops off with 
height as a potential field we thus obtain 

Q0) cx J 2 (z) oc B 2 (z) cx [exp (-fc z)f 

oce- 2feoZ , (1) 
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Fig. 8. — Similar to Fig. [S] here for loop 2 with many reversals in gas flow direction 



and so the estimated heating scale height is sh — 
l/2fco ~ 5.0 Mm, consistent with what we actually mea- 
sure. 

In terms of th e half loop length L of a ssumed semicir- 
cular loops (cf. lAschwanden et afll200H) . our estimated 
heating scale height is 



where the wave number fco and the half loop length L 
have been related to the distance between polarities d 
through fcn = ir/d and L = nd/4. This is consistent with 
lAschwanden et al.l (|2001^ . who find s H /L = 0.2 ± 0.1. 

3.4. Loop stratification and dynamics 

Loops encountered in this model are generally not in 
equilibrium and a significant fraction of the visible loops 
identified in the TRACE 171 filter are not aligned with 
the magnetic field. A possible scenario for such a loop 
to develop is when one of two connected polarities has 
a velocity gradient in the direction connecting the two 
polarities. This implies that one polarity is drawn out 
along a line perpendicular to the line initially connecting 
the two polarities, which will create a shear in the field 
connecting the two. In general the maximum of the shear 



will not follow a particular fieldline connecting the two 
polarities. If the shear happens on a short timescale the 
maximum shear locations will enter the narrow TRACE 
filters at the same time, showing a bright region that is 
not connected by a single loop bundle. 

Three loops are traced in Fig. 21 where loop 1 and 
loop 2 are selected for their brightness in the emulated 
TRACE 171 filter (Fig. H top) while loop 3 is selected 
for its brightness in the emulated TRACE 195 filter (Fig. 
H middle). 

Loop 1 can be identified through its whole length in 
the emulated TRACE 171 filter. It is an example of a 
steady relatively cold loop and seems to be in a cooling 
phase. It shows low steady velocities down both of its 
legs of at most a few tens of km s _1 . The energy balance 
and the gas parameters of the loop are shown in Figures 
[S] and respectively. The heating is growing smoothly 
towards the photosphere, except for a small bump on the 
right leg, which could be an impulsive heating event. The 
convective energy term advects energy from the coronal 
part of the loop to the denser lower atmosphere, were it 
is deposited. The Spitzer conductivity shows the charac- 
teristic pattern of moving energy from the hot corona and 
pumping it down into the chromosphere. The radiative 
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Fig. 9. — Similar to Fig. [S] here for loop 3 
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Fig. 10. — Similar to Fig. |5] here for loop 3 with gas flows down both legs from a position just below the loop top in the left leg 



cooling is insignificant and shown here without quench- 
ing, making it rise sharply when the density increases. 
This loop is very quiescent except for the bottom of the 
right leg, which shows large fluctuations in the convective 
term, due mainly to the large densities, and not so much 
large velocities. The temperature, density and pressure 
are very close to constant in most of the corona. There 
is a slight decrease in temperature, which is large enough 
to have an effect on the emission in TRACE 171 filter. 

Loop 2 is shown in Figures [7\ and |SJ This loop is an ex- 
ample of a dynamic loop that, in spite of the intermittent 
velocities present in the loop, is still visible throughout 
most of its length. This loop has more uneven heating 
than loop 1. The velocities present in this loop are the 
effect of the heating and velocities perpendicular to the 
loop. This loop is much longer than loop 1, even though 
this is not obvious from Fig. 01 (middle), the top view 
shows that this loop is much longer than loop 1. Pro- 
jection effects such as this one arc quite common, and 
lead to confusing and at times misleading conclusions 
about a particular loop. The special geometry of the 
loop, with a flat loop top, makes the gradients in veloc- 
ity in Fig. [S]seem strong because the velocity is plotted 
against height. In spite of the gradients being smaller 
than what they seem from Fig. [7| there is still a signifi- 



cant energy reorganisation in this loop. The convective 
energy term is dominating everything else in the corona 
produced by velocities of up to 30 km s _1 in the right 
leg. The left leg has a distance over which the velocities 
are on the order of 1 km s _1 . The temperature along 
the loop is close to constant through the whole corona, 
while the pressure and density are larger in the left leg, 
with a noticeable increase in density at the loop top. It 
is because the loop top is flat that it is possible for this 
amount of mass to have accumulated at the loop top. It 
could have happened as the result of powerful heating 
at both legs, which has now turned off. The evaporated 
gas would then accumulate at the loop top, and now the 
loop is cooling down. The low pressure of the right leg is 
at this point in time not able to support the dense gas at 
the loop top, which is now falling down the right leg. At 
the same time there is still material coming up the right 
leg, and a large amount of convective energy is deposited 
where the two flows meet. The left leg also has material 
coming down from the top, but this is slowed down high 
in the corona, with a corresponding effect on the convec- 
tive energy term. Below this point gas seems to be free 
falling due to the very low magnetic heating in the left 
leg. 

Loop 3 is somewhat hotter than loop 1 and 2, and 
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shows up most prominently in TRACE 195. Figures 
and llOl show that the heating in the right leg is more than 
a factor of ten higher than in loop 1 and 2, and further- 
more is much larger in the right leg, compared to the left. 
From Fig. 0] one can see that the right leg is rooted in 
the main polarity, with a magnetic field strength roughly 
ten times higher than the left leg. Because of the high 
pressure in the right leg, the transition region is also at 
different heights in the two legs. In the left leg, the tran- 
sition region is at a height of 3.5 Mm while in the right 
leg it is at 7 Mm. As already mentioned the heating is in 
general dependent on the magnetic field strength and the 
local height of the transition region. Because of the large 
differences in magnetic field strength and local height of 
the transition region, the heating is much stronger in the 
right leg than in left. In spite of the powerful heating 
in the right leg, the radiative cooling close to balances 
the powerful heating. The pressure in the right leg is 
building, and is now higher than at the top of the loop. 
Nevertheless, the momentum of the falling gas is still 
large enough that the increasing pressure gradient is not 
yet able to break the fall. 

4. TIME EVOLUTION 

Loops observed with TRACE and SOHO/EIT can be 
long lived. Some show life-times much longer than dy- 
namical, cooling and conductive time scales. In following 
a loop in time we have attempted to follow an emitting 
plasma-parcel which initially is caught in a flux bundle. 
The plasma parcel is then followed to the next snap shot, 
and here the magnetic field is used as a tracer for the rest 
of the loop. In a 3D simulation like the one discussed 
here, following a plasma-parcel from one snapshot to the 
next is theoretically simple. But following a plasma par- 
cel over a long time may require sub-resolution precision, 
because errors in position accumulate from snapshot to 
snapshot. Unless the plasma parameters change abruptly 
there is no way to estimate if an error has been made 
when following a plasma parcel over long times. In spite 
of these problems we believe we have managed to follow 
the time evolution of loop 1 with good precision for more 
than 20 minutes. The emission in TRACE 171 is shown 
in Fig. ^2 m units of DN s _1 pixel -1 Mm -1 (making 
the actual photon count rate depend on the thickness of 
the loop as well as on the background and foreground). 
One resolution element, emitting at the level of most of 
the upper part of the loop, makes up 10 % of the typ- 
ical level of emission from a loop seen in Fig. The 
loops in this simulation are in general thicker than in the 
solar corona, because of the limited resolution in the pho- 
tosphere, which sets the minimum driving scale. After 
the expansion of the magnetic field through the chromo- 
sphere and transition region, this makes the loops in the 
simulated corona thicker than on the Sun. In this simu- 
lation, loops generally are on the order of 6-8 grid cells 
thick or 3 Mm. In spite of the very dynamic nature of 
the loop 1, it manages to stay bright along most of its 
length during the whole 20 min time period. The dark 
left leg of loop 1 is caused by a too high temperature, 
moving that leg outside the temperature response func- 
tion of the TRACE 171 filter. The reason that the leg is 
still apparent in Fig. 01 (top) is the large amount of other 
loops originating in the same area. Fig. II II clearly shows 
that it is possible to have long lived, non static loops in a 
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Fig. 11. — Emission in the TRACE 171 filter as function of time 
for loop 1. The dashed line shows the time of Fig. 1^1 



corona produced by a highly time dependent DC heating 
mechanism. 

5. DIFFERENTIAL EMISSION MEASURE 

DEM can be constructed from observations of a num- 
ber of lines, essentially giving a measure for the amount 
of emitting gas at a certain temperature. DEM can 
also be constructed by comparing emission in a num- 
ber of narrow filters, such as the ones used by EIT 
and TRACE. Constructing DEM curves on the basis 
of narrow band observations is however an ill posed 
problem, and whether any trustworthy information can 
be gained from suc h exercises has been the subject 
of intense discussionllSchmelz et al.ll200H iMartens et alJ 
12003 lA^chwandenll2002tr Creating DEM curves from a 
3D simulation is more straight forward. The DEM can 
be defined as 

DEM = n e ^, (3) 

which, when defining a volume emission measure (EMy) 
as 



EMy(T) 



nldV . 



(4) 



makes it possible to find the DEM from the simulation 
simply as 

DEM = -r-T^— (5) 
A AT w 

where A is the area of a grid cell. When constructing 
DEM curves AT is often not constant, but instead con- 
stant logarithmic intervals are used. In this case we use 
temperature intervals of dLog 10 T = 0.1. The intervals 
are furthermore over sampled in the sense that inter- 
vals overlap in order to be able to construct reasonably 
smooth curves. 

DEM curves for 5 selected points along loop 1, shown 
as crosses in Fig. 0] (top), are shown in Fig. El It is very 
difficult to identify loop 1 from these curves, and it is also 
clear that especially for points 2 and 3, the typical tem- 
perature visible in TRACE 171 is not the temperature 
where most of the gas is. Point 5 shows especially many 
isolated bumps, while point 3, even though it seems to 
be in the middle of a bright loop, show a clear bump 
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Fig. 12. — DEM curves for the five points along loop 1 shown 
in Fig. [U The points are numbered from left to right with num- 
bers 1 to 5. Curves 2-5 are shifted up by 1-4 orders of magnitude 
respectively 



at low temperatures. These DEM curves are representa- 
tive of most of the ones we have produced and one can 
only seldomly find clear signatures of loops, because of 
the "pollution" of loops seen in images like Fig. 0] (top). 
In spite of these problems, it seems that in this model, 
loops typically have a DEM of 10 21 - 10 22 cm~ 5 KT 1 near 
1 MK. This level of emission measure is also what for 
instance iSchmel z et alJ (l2003f) finds. In the same work 
(|Schmelz et al.ll2003|) argue that because the maximum 
for their DEM curves along the loop is not at the same 
temperature one has to be very careful in employing this 
method. We agree on this point, but not on the conclu- 
sion that the non constant maximum of the DEM curve 
is evidence for a temperature gradient along the loop. 
The conclusion would hold true if the loop were per- 
fectly isolated in temperature. This is very seldomly the 
case. By visual inspection, loop 1 seems to be reason- 
ably isolated in Fig. 0] (top) , but we have seen that at 
the left leg this is not so, even though along some of its 
length it is true. Figure shows the temperature of the 
gas and the contribution to the emission in TRACE 171 
along the projected axis in Fig. Q] (top) of point 4. It is 
also clear that the sharp peak in emission is located in a 
temperature dip, and that even though there are quite a 
few areas at a temperature which should show up in the 
TRACE 171 filter, these are not important, because the 
density is too low. The DEM is able to show the range 
of temperatures involved but in spite of the sharp peak 
in emission in the TRACE 171 filter, it is not clearly 
noticeable in the DEM curve. It therefore seems that a 
simplistic approach to DEM curves is a problematic way 
of trying to identify loops. 

The problem would become more severe when also the 
resolution is worse than the typical width of a loop, since 
this would introduce not only contamination along the 
line of sight, but also in the 2D viewing plane. 

6. DISCUSSION AND CONCLUSIONS 

The three loops discussed here are representative of a 
large number of loops we have examined and illustrate 
that there are no static loops. Loops are continuously 
changing, with large differences in heating, cooling and 
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Fig. 13. — Emission contribution (solid) and temperature 
(dashed) along the projected axis at point 4 of loop 1 shown in 
Fig. H (top) 



cross sections. It is therefore futile to fit scaling rela- 
tions to individual loops, since the scaling properties will 
change with time and between loops, making any scaling 
relations inexact. We have seen examples of loops in the 
emulated TRACE pass bands that do not follow field- 
lines, because the field is sheared and heating is localised 
in thin structures across fieldlines. The high thermal 
conductivity adjusts the temperature gradient until con- 
duction balances the heating rate, and makes the effect 
hard to observe except when projection effects add up 
contributions from the slightly brighter sheared heating 
region. 

When looking at large ensembles of loops in an active 
region the heating scales roughly with the magnetic field 
strength squared and decreases exponentially from the 
local height of the transition region. The scale height 
is essentially proportional to the mean separation of the 
main polarities in the active region. 

Velocities in loops are assumed to be along fieldlines, 
and that is what we find in the region just above the 
transition region. Above roughly 10 Mm the velocities 
along the loops relative to the total velocity amplitudes 
are no longer close to unity, but span the whole interval 
between zero and one. The velocity amplitudes are on 
average a few tens of km s _1 , but velocities as high as 
400 km s^ 1 are present. 

Density and pressure are less extreme in these 
loops than what has b een reported for long loops by 
iWinebarger et all ()2003[K and are here only a few times 
denser than the surroundings. It is, however, clear that 
the density is approximately constant along loops. The 
loops therefore show no decrease in emission measure 
with height. This simulation only spans 37 Mm in height, 
slightly less than a pressure scale height for a 1 million 
degree gas, and therefore we cannot say if higher loops 
would show a decrease in density with height but since 
the heating scale height also tends to increase with ac- 
tive region size, we expect loops to be able to keep ap- 
proximately constant density, or at least have densities 
decreasing slower with height than for an atmosphere in 
hydrostatic equilibrium. 

Both heating and plasma flows are time dependent for 
loop 1, even though it has an almost constant emission 
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throughout the time span we followed it. The coronal 
loops analysed here may be assumed to be generic, but 
are most likely not representative of many of the loops 
observed by TRACE and analysed in detail elsewhere 
since the latter are biased towards being bright, isolated 
and long. The DEM curves measured for the bright 
TRACE 171 loop do not easily reveal a loop, and as 
already concluded bv ISchmelz et alJ l)2003j) . differential 
emission measures must be treated carefully if they are 
to be used. 

The magnetic field in this work is only disturbed 
slightly from a potential field, so large scale helicity is not 
present. Large scale helicity is often present in observed 
active regions on the Sun and may be caused by large 
scale velocity fields, or may exist already in the flux emer- 
gence phase. Including disturbances able to create sys- 
tematic helicity would introduce unknown or arbitrary 
variables in contrast to the ab initio approach of this 
work, but we do expect that if a number of short loops 
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in quiescent active regions are observed in the TRACE 
or EIT 171 band similar loops as modelled here would 
be present. 

Analysis of this simulation indicates that often-made 
assumptions, such as those of hydrostatic equilibrium, a 
time independent and symmetric heating function, and 
constant cross sections along loops are highly question- 
able and cannot provide a basis for simplified models of 
loops. 
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